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Abstract 

We study the process of inhomogeneous tachyon condensation in an inter- 
secting Dl- and anti-Dl-brane system using an effective tachyon DBI action. 
By switching to the Hamiltonian formahsm, we numerically solve for the dy- 
namical evolution of the system at a small intersection angle. We find that 
the decay proceeds indefinitely and resembles the action of two zippers mov- 
ing away from the intersection point at the speed of light, zipping the branes 
together and leaving inhomogeneous tachyon matter behind. We also discuss 
the range of validity of our analysis and discuss the relation of the Dl-anti-Dl 
description of the system to one in terms of an intersecting Dl-Dl-brane pair. 
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1 Introduction 



The decay of branes via tachyon condensation is one of the few tractable string theory 
systems with true nontrivial, nonperturbative dynamical evolution. The process by 
which branes decay or annihilate is of fundamental relevancy and reveals important 
connections between open and closed string sectors. While much recent progress has 
been made in understanding open string tachyon physics, the system is still not well 
understood and many open questions remain. 

In particular, the annihilation of parallel D-branes and D-branes is sufficiently 
simple to allow both a limited worldsheet analysis [1] as well as an effective descrip- 
tion in terms of a single homogeneous tachyon field [2,3]. The inhomogeneous case, 
however, is both more interesting and more complicated. For example, chiral sym- 
metry breaking in the Sakai-Sugimoto model is an example of a physically relevant 
localized tachyon decay [4,5]. Much of the effort so far has been limited to study- 
ing time-independent inhomogeneous soliton solutions and marginal deformations [6] 
rather than on the dynamical evolution. 

Perhaps the simplest inhomogeneous system with a localized tachyon is where two 
D-branes intersect at an angle 9. The tachyonic mode is localized at the intersection, 
and as it condenses, the two branes reconnect. This recombination process plays 
a central role in realizing the Higgs mechanism in Standard Model on intersecting 
D-branes [7] and cosmological models of brane inflation [8]. In addition, this type 
of localized condensation may serve as a toy model for the reconnection transition 
found, for example, in the Sakai-Sugimoto model, in non-trivial curved backgrounds. 

For any nonzero value of the angle 9, as mentioned above, the open strings stretch- 
ing between the branes have tachyonic modes in the low-lying mass spectrum. As is 
well-known from worldsheet calculations, the negative mass squared of the tachyon in 
the NS (Neveu-Schwarz) sector increase linearly with 16*1 [9]. It is important to notice 
that when the angle 9 becomes large, these are no longer the lowest mass excited 
states of the system. In particular, when the intersection angle 9 passes through | 
(perpendicular branes) a new tower of states becomes increasingly light and constitute 
the low-lying mass spectrum of the system [9-17]. Indeed, for the most extreme case 
~ vr, it is more natural to view the system as a D-brane and a D-brane intersecting 
at an angle (p = -k — 9. 

Ideally, one would like to study the intersecting D-D system directly on the world- 
sheet by calculating the string scattering amplitudes in the rolling tachyon back- 
ground [2, 18]. However, this requires turning on both the inhomogeneous tachyon 
deformation and a transverse scalar at the same time, a task which is currently out of 
reach (see [1,19,20] for work toward this direction in different set-ups). We therefore 
resort to describing the inhomogeneous tachyon condensation from an effective fleld 
theory point of view. 

For economic reasons we will consider here intersecting Dl-branes. The general- 
ization to other Dp-branes and multiple angles is straightforward. 

In both flat and curved space, there are known constructions for coincident D- 
D effective actions [3,6] and the consequent work on understanding their homogeneous 
decay [21]. However, when D-branes are even slightly off from being parallel to each 



2 



other, little is known about their decay process in real time. This is mostly because 
initially the branes meet at an angle at a single point in spacetime making the tachyon 
condensation process highly inhomogeneous. 

To derive an effective action for a Dl-Dl system, we will start from a known action 
for a coincident D2-D2-brane pair with appropriate gauge fields and tachyons turned 
on. Then, after T-duality, we will be led to an action describing a Dl-Dl-brane pair 
which is initially intersecting at an angle ip in flat spacetime background. An analysis 
of the spectrum of small fluctuations around the false vacuum will give a tachyon 
mass which, by construction, is exact for small y?. 

Our main objective is to solve for the explicit temporal evolution, one which 
we will have to approach numerically. The inhomogeneous decay is described by 
two coupled two-dimensional fields whose equations of motion are not amenable to 
analytic methods. While we can not entirely trust the quantitative results, because 
we are solving an effective theory and in addition doing so only approximately, our 
solutions will hopefully correctly capture the qualitative features of the decay. 

The branes decay in roughly two steps. First, starting from the intersection point, 
the separated Dl and Dl are attracted to each other in a manner that resembles being 
zipped together with the zipper traveling asymptotically at the speed of light. Once 
the branes are parallel, the tachyon begins to roll with constant velocity toward its 
vacuum, and, just as in the homogenous case, the dynamics of the branes can then 
be described in terms of a pressureless, non-interacting tachyon matter. 

Our results qualitatively match the expected evolution of straight intersecting 
branes reconnecting into curved separating branes, as has been seen in previous anal- 
yses in terms of Dl pairs [15,17], except that in this Dl-Dl description the separating 
branes are connected by a parallel Dl-Dl pair with a rolling tachyon. We will ar- 
gue, using a slightly modified example, that from the open string point of view, this 
parallel Dl-Dl pair is equivalent to the vacuum, and so that the final states of the 
Dl-Dl and Dl-Dl systems in fact differ only by a change of variables. 

This paper is organized as follows. In Section [2] we will review previous work on 
tachyon condensation in intersecting D-brane systems, and in Section [3] we will turn 
to deriving an effective action for Dl-Dl-branes. In addition, we will compute the 
false vacuum modes and the tachyon mass. Section H] contains some steps necessary 
for numerical evaluation of the dynamics including a Legendre transformation to the 
Hamiltonian formalism and a discussion of the boundary conditions. In Section [5] we 
will present and describe the results of the numerical computations. We will then 
explain in Section [HI the relationship between describing the decay as that of two 
D-branes rather than as a D-D pair. Finally, Section [7] comprises our summary and 
ideas for possible extensions of this work. 

2 Review of previous analysis 

Let us first recall some basic facts from the worldsheet. For simplicity, let us start 
with two parallel Dl-branes in Type IIB superstring theory. Then consider dialing the 
intersection angle 9 to some finite value. The open strings which connect the different 
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Dl-branes will then get confined about the intersection point because the tensions of 
the strings tend to minimize. By analyzing the NS sector of the energy spectrum of 
these localized strings one finds that the mass squared of the modes behaves linearly 
with the angle [9], 

1\ 1^1 



'^ws=(N-- j-^ , iV = 0,l,2,... . (1) 

In the rest of the paper we shall only focus on positive angles, ^ > and (p = tx—9 > 0, 
and will henceforth drop the absolute value symbols. We will be interested in the 
lowest, tachyonic excitation = 0. Notice that the tachyonic mode exists at all 
(nonzero) angles 6 and thus signals the instability of the configuration. 

The part of the spectrum in ([T]) is reproduced to order 0{9) by the spectrum of 
fluctuations around the intersecting Dl-Dl-brane pair background in a non-Abelian 
Yang-Mills (YM) theory [15,17]. An intrinsic feature of the effective field theory 
approach is that the mass of the lowest mode behaves as 

, tan ) ^ ^ 

^2 ^ ^ ^ (2) 

rather than linearly as in ([T]). From ([2]) it is evident that the spectra only match 
at small angles, but the situation becomes increasingly worse at larger angles. This 
reflects the fact that the YM description is only viable for small angles. The negative 
mass squared of the tachyon means that the eigenfunction blows up exponentially in 
time [15,17], 

tan ( ^ J 

T ~ e-'^*e-^^" . (3) 

From (I3j) we also see that the tachyon fluctuation modes are localized around the 
intersection point x = 0, agreeing with the worldsheet. In [15] it was further argued, 
that the geometric realization of the tachyon condensation is a D-brane recombination 
process. This was explicitly shown by diagonalizing the fluctuations through a local 
gauge transformation of the brane-coordinates. Here, however, we shall postpone 
discussing this phenomenon to Section [6l 

One can also use the spectrum ([T]) to check the a'^F^ and higher order a' correc- 
tions in the expansion of the non-Abelian Dirac-Born-Infeld (DEI) action and find 
agreement [11,16,17]. However, equipped even with the full DBI action, thus far one 
has not been able to reproduce the mass spectrum ([T]) exactlyQ 

In the current paper, we are mostly interested in the behavior near ^ ~ vr. Given 
the shortcomings of matching with the mass spectrum ([T]) it is better to look for 
alternative routes. If one dials the intersection angle all the way to ^ — > tt, one finally 
arrives to a parallel Dl-Dl configuration, which has been investigated in numerous 
articles. We will therefore find it promising to begin with a well-motivated tachyon 
DBI action for a Dl-Dl-brane pair, which, by construction, will be valid at small 



^In [11] it was attempted to fix the discrepancy by considering a symmetric trace prescription 
of [22] in the non-Abehan DBI action. Though this procedure did fix the behavior for the lowest 
excitation, the spacing of the mass spectrum turned out to be incorrect. 



4 



angles (p = tt — 6 ^ 0, in contrast to the usual effective action for the Dl-Dl-brane 
pair. Derivation of the Dl-Dl action shall follow below. 



3 Derivation of the Dl-Dl action 

We begin by deriving the effective action for the intersecting Dl-Dl system. This 
can be acquired from the action for a coincident D2-D2-brane pair with equal and 
opposite magnetic fields. We then perform a T-duality along one of the directions 
longitudinal to the D2 and D2 which yields the desired Dl-Dl action. 

The action for a pair of coincident D2-D2 is the usual tachyon DBI type action [6], 



Sd2 = -/i2 I d'x Vi\T\) [e-^'^y-det (Xl^) + (1 ^ 2) j , (4) 



where 

4? = r/,fe + 9„Z»a,Z«+2WF« + ^(Z},r(Ar)* + (Z}„r)*Ar) (5) 

f!S = 0.4^ -d,A^^ (6) 
D^T = 9,r-z(A«-Af)r. (7) 

By choosing to work in unitary gauge, we can set Im r = and will denote Re r = T. 
With a view towards an upcoming T-duality in the y-direction, let us choose an ansatz 
where T and Ay depend on t and x but not y. We further set A^'^ = —A^^ = A. In 
addition, we fix the dilaton to be constant e*^''' = gs and set the other scalars Z^^^ and 
the other components of the gauge field to zero. The dynamical variables are then: 



(9) 
(10) 

(11) 
(12) 



Fty 


= dtAy = A 


J- xy 


= dxAy = A 


DtT 


= doT = f 


DxT 


= d.^T = T' 


DyT 


= -2iAT 



With this ansatz, D2 and D2 actions are identical and combine into 

Sd2 = -2^ [ d^'x y(T) V- det (r/,, + 27ia'Fab + 27ra'DaTD,T) . 

gs J 



(13) 



Evaluating the determinant and integrating over the circle in the ^/-direction of radius 
R yields 



Sd2 = -2(27ri?)^ j d\V{T) y(l + 2W (-T^ + t")) (l + SWA^T^) 

+ {2Txa'f (-A^ + A'') - (2^)=^ {Ar - A'f^ . (14) 
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We now T-dualize in y to find the action for intersecting Dl-Dl. Under tlie 
T-duality, R ^ gs ^ "^5'*' "^v ~^ 2™^^' "^^^^ yields tlie action: 

So, = -2^ [ d?x V{T) (15) 

gs J 



X W (^1 - 27ra'T2 + 2Txa'T'^^ {l + - v"^ + V'^ - 27ra'{y'f - yT] 
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wliere fxi = 27r-\/a'/i2- The scalar y represents half the distance between the Dl and 
Dl, so the initial angle between the branes is then given by 

tan(^/2 = - . (16) 

X 



We can rescale the coordinates t and x and the field y by {t,x,y) \/2na'{t, x, y), 
so we have an action in terms of dimensionless quantities: 



= -J^ j d^x V{T)^{l - t2 + (1 + 4T2?/2) _ ^2 + y,2 _ ^yif _ yj^iy ^ 

(17) 

where we define the normalization Af = 2(27ra')^ = y. Motivated by [23], we take 
the tachyon potential to be 

= ^Jm ■ 

where (3 = -^/tt for superstrings. The vacua for this potential are at T = ±00, which 
implies that to reach them would require an infinite amount of time. 

3.1 EM tensor for tachyon matter 

For later purposes let us record the energy-momentum tensor for the action fll7p . Let 
us generalize the above action (fTTIl to curved space as follows 

^curved = j dtdx^ V {T) ( 1 + g''\dayd,y + daTd,T) 



+ATY{l+g'^'daTd,T) + {g''''g'''-g'''g"')dayd,yd,Td,T\ . (19) 



The energy-momentum tensor is then extracted as 

2 (^i^curvcd 



'-ab 



-g Sg 



ab 



(20) 



MV{T) 

VabJL' + 



(^1 - r2 + T'2j (1 + 4r2y2) -y^ + y>^ - {y'f - yVf 

X ( daydi,y + daTd,T + AT^daTd^T + dayd,y{dTf 

+daTd,T{dyf - d(^aydb)Tdy ■ 9t) . (21) 

6 



Explicitly, 



iiv' + TT' + AT^v'^TT' 



(^1 - T2 + 2-/2 j + 4T2y2) - y2 + ^/2 _ (^^/^^ _ ^2-/)2 



_1 + y2 , _ 42-2 2/^ _ 2-2N 



(^1 - T2 + 2-/2^ (1 + 4T2|/2) - y2 + ^/2 _ ^yff _ yrpiy 



3.2 Mode analysis 

As in the previous studies of D-D systems [15, 17], we can analyze the spectrum of 
fluctuations of the action ( ITTI) around the initial state to identify the tachyonic mode. 
We begin at the maximum of the tachyon potential T = and will also allow for the 
possibility of small x-dependent corrections to the relative position of the Dl and Dl, 
so y{x) = xtan{ip/2) + Sy{x), where 6y <C 1. The action (|T7I) expanded to second 
order in T and 6y is then 

S ~ Ary"rfx^-^ + ^ + ^(4x2tan(y./2)2-/?2(i + tan((^/2)2))r2 

+ t8.n{^/2)6y~{6yf + iSy)'^ ■ (25) 

The position fluctuations 6y are free and decouple, so at this order the position of 
the branes is uncorrected. 

The equation of motion for the tachyon is 

-f + T" = (4x2tan((p/2)2-/52(i+tan(v^/2)2))r . (26) 

Now, decompose T into modes of definite frequency Qn as 

oo 

T{t,x) = Y,Cne~'''"'T4x) . (27) 

n=0 

Normalizability of the modes implies a discrete spectrum, and therefore n is a non- 
negative integer. The equation for the tachyon fl26p now becomes just the Schrodinger 
equation for a harmonic oscillator with mass m and frequency u, where mu = 
2tan((/?/2) and 2mEn = /5^(1 + tan(y9/2)2) + Qf^. Imposing the boundary condition 
T„(oo) = 0, the solution for the modes is therefore 



T„(x) = Hr, (xv/2tan(<^/2)) e"^' ^^^^'^/^^ , (28) 
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where Hn is the Hermite polynomial of order n. These modes are localized, as ex- 
pected, near the intersection point at a; = 0. 
The frequencies are given by 

nl = 4tan((^/2) + ^ - (l + tan{^/2f) . (29) 

Focusing on the lowest mode n = and plugging in (3"^ = tc, we see that Qq < for 
all values of ip, giving an exponentially-growing tachyonic mode. For small (f, 

nl^ip-7i (30) 

which matches with the worldsheet calculation ([T]) of the tachyon massH However, 
as with the D-D calculations in [15, 17], for larger angles the mass computed via the 
effective action increasingly deviates from the exact result. 



4 Set-up 

Having arrived at the effective action (fT7|) . we now need to derive the Hamiltonian 
equations of motion. In addition, we will also discuss the appropriate boundary 
conditions, particularly the modifications required for implementing numerical com- 
putations. 



4.1 Hamiltonian formalism 

We will perform a Legendre transformation on the action ffT7|) and work in the Hamil- 
tonian formalism [24-26]. Of course, the Lagrange equations of motion derived from 
the action ( ITTI) are in principle equivalent, but for implementing a numerical solu- 
tion, solving four coupled first-order differential equations was found to be easier than 
two coupled second-order differential equations. In addition, using the Hamiltonian 
framework facilitates a description of the tachyon vacuum, V(T = ±oo) = 0. 
The canonical momenta are given by 

dC -^ViT) {t (1 + Ay'T') + y' [yV - y'f) } 

Hy — — = (31) 



dT 



(^1 - r2 + T'2j (1 + 4T2y2) - ?/2 + y'2 _ (^y'f - yVf 



dC UV{T){y + T'(yr-y'f)] 

n, = , I V yj (32) 

y ' (^1 - T2 + 2-/2 j (1 + 4y2^2) _y2j^ y,2 _ (^y,f _ yj^^y 



^Recall that we are measuring dimcnsionful quantities in units of 2Tra' . 



8 



The Hamiltonian density is then 

H = Urf + Uyy-C (33) 

= ^n|(l + T'') + + 4y^T' + y'') + 2y'T'n,nr 

1 

+M''V{Tf {y'^ + (1 + r'')(l + Vt^)) j (34) 



from which we derive the equations of motion: 
. _ IiT{l + T'^)+y'T'ny 



(35) 



n 

-myr - M^virf {4 (i + r^) y^r + ^ {y'^ + (i + r^) (i + 4tV))} 



( ^ir + y-n^n, + ArV(T)^r (i + AT^y^) ^ ^^^^ 

. _ ny{l + Ay^T^ + y'^)+y'T'nT 



-4n2,T2|/ - AU^V{Tf (1 + T'2) T^y 



+4. 1 a^™%±^™: , , ,33) 



7^ 

n 

~n 

4.2 Boundary conditions 

For describing the dynamical evolution, we begin at t = with a straight, static 
Dl and Dl intersecting at an angle ip. This translates to ?/(0,x) = xtan((/?/2) and 
y(0,x) = 0. In order to initiate the tachyon rolling, we start it at rest, T(0,x) = 0, 
but slightly displaced from the maximum at T = near x = 0. We do not consider a 
homogenous tachyon, both because the decay is localized near the intersection and, 
more importantly, as we will describe, for performing the numerical calculations the 
tachyon must be zero at the spatial cutoff Xmax- The initial tachyon profile 



2 2 



T(0, x) = T, (^e-^ - e"^— j (39) 

is somewhat arbitrary, but we have checked that the subsequent evolution is largely 
insensitive to its shape. At the level of small perturbations at least, we can decompose 
the initial profile into modes, and while the tachyon mode grows exponentially, the 
massive modes will just oscillate, contributing some small wiggles. Changing the 
overall constant simply adjusts the timescale of the decay, as we will discuss in the 
next section. 

Since we will be solving the Hamiltonian equations of motion, we need to translate 
our initial conditions for T and y into conditions for 11^ and 11^. From (l3Tl) and (j32l) 
we find that Urit = 0, x) = and Uy^t = 0, x) = 0. 
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In addition to initial conditions, we need to impose conditions at the spatial 
boundaries. Ideally, we would like to require that at spatial infinities the system is in 
the false vacuum for all times. However, for numerical calculations, we employ a finite 
spatial cutoff, x G [— a^max? ^^max]- At these boundaries we need to fix the fields, and 
they must satisfy the equations of motion there. We therefore fix the finite boundary 
to be in the false vacuum, i.e., T(t, ±Xmax) = and y{t, ±Xmax) = ±a^max tan{ip/2). Of 
course, we must be sure to choose Xmax large enough that the spatial cutoff does not 
affect the dynamics. Since boundary contributions can only propagate at the speed 
of light, the region outside the future lightcone of the boundary should be insensitive 
to the cutoff. We will be focussed on the relevant physics near the intersection point, 
so if we consider only times t < x^ax we will be free of spurious boundary effects. As 
a check, we have varied Xmax and shown that for sufficiently large values, our results 
are independent of the choice. 



5 Analysis 

We solve the Hamiltonian equations ( l35l) . ( 136|) . ( 1371) . and ( l38l) using Mathematica. 
For the results presented below, we choose the following representative parameters: 

= 1 (40) 

Xmax = 30 (42) 
T, = 10"^ . (43) 

We find that the numerical integration is only trustworthy up to some value of t. 
After this point numerical errors begin to grow large; for example, the stress tensor 
is no longer even approximately conserved. For these parameters, we found that the 
numerics broke down around t > 10. 

The solution for the tachyon as a function of t and x is shown in Fig. [H Initially, 
the system remains close to the false vacuum. But then, near x = 0, the tachyon 
begins to grow, and the region of condensation spreads outward roughly at the speed 
of light. The corresponding evolution of y is shown in Fig. [21 For clarity, because 
y is antisymmetric, y{t,x) = —y{t, —x), we only plot positive values of x. Initially, 
the system again stays near the false vacuum and the branes are straight. Once the 
decay process begins, for x <t the branes are coincident with some small oscillations 
around y = 0. 

We can more easily illustrate certain features of the evolution by considering 
cross-sections of Figs. [U and O In Fig. [3] we plot T(t,0). There are three clearly 
distinguishable parts of the curve; at early time, for t < 3, the system remains close 
to the false vacuum T ^ 0, then there is a brief transitional region where T begins 
to condense, rolling down the potential, and then for t > 4, the tachyon rolls at a 
constant velocity T ^ 1 toward the vacuum at T = oo. 

In addition, we plot in Fig. H] the profile of y{x) for various fixed times, again only 
for X > 0. For t = 4, the branes are just beginning to pull toward each other. At 
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Figure 1: The tachyon as a function of t and x. 





Figure 4: y as a function of x, for t = 2 (black), 4 (dotted red), 6 (dot-dashed blue), 
and 8 (dashed green). 



zipper 




Figure 5: The position of the positive- velocity zipper z as a function of t. 

later times, the profile assumes a waterfall shape, again with three regions. Far from 
the intersection the branes retain their initial configuration. But, at around x ~ t, a 
point is reached where y — > rapidly, and then, in an increasingly large region around 
X = 0, the branes have become coincident except for some small fluctuations. The 
process resembles the action of two zippers moving away from the intersection point, 
zipping the branes together. 

We will define z{t) to be the position of the right-moving zipper and by symmetry 
the left-moving zipper is at —z{t). More specifically, for a given time, z equals the 
largest zero of y{x). Fig. [3] presents a graph of z{t). For t < 4, z = while at late 
time z ^ 1 with a sharp transition in between. 

We can further illustrate the physical properties of the evolution with plots of 
the components of the stress tensor. Fig. [6] shows the energy density Too. We have 
normalized the false vacuum energy so that Tqq = 1. The motion of the zippers 
matches roughly with the large spikes in the energy density due to the large kinetic and 
gradient energy as the branes come together very rapidly. The position of the zippers 
is also clearly visible in Fig. [7] which plots the momentum density Tqi; the zipper 
moving to the left is represented by the spike and the one going to the right by the 
trough. The pressure Tn is shown in Fig. [HI The false vacuum has negative pressure, 
and again there are spikes corresponding to the worldlines of the zippers. Once the 
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Figure 6: The energy density Tqq as a function of t and x. Notice the large spikes 
corresponding along the worldlines of the zippers. 

tachyon has begun to condense, Tn 0. In addition, we checked numerically that 
Tf^i, was conserved to very good accuracy by the evolution. 

We will first focus on the transitional region between the false vacuum and the 
rolling region, where the decay first begins. We can identify three processes which 
characterize this region, and we can see that they are all related. Intuitively, the 
tachyon is massive where the branes are separated by more than the string length, so 
it cannot really begin rolling until the branes have been zipped together. However, 
it is the displacement of the tachyon from zero which gives a potential to y and 
pulls the branes together. Consequently, both of these processes must occur together. 
The timescale for the onset of the decay is controlled by the initial displacement 
parameter T^; the smaller is, the longer it takes for the decay to begin. For our 
choice of = 10~^, the tachyon starts rolling near the intersection at around t ~ 3. 
At later times T{t, x) starts to grow when t ^ x. Similarly, the zipper begins moving 
at t ~ 4. As the zipper's velocity approaches the speed of light, a given position on 
the brane begins to deviate from the false vacuum at around t ~ x. 

The late time behavior, once the zipper has passed, very much resembles the 
decay of parallel Dl-Dl-branes but is inhomogeneous. We actually do not find that 
y = past the zipper but rather that y overshoots and oscillates around zero. As T 
increases, the masses of these oscillations grow, so their frequencies and amplitudes 
diminish. However, we should not really trust our effective action to describe features 
such as these wiggles which are small compared with the string length. 

In the approximation that at late times y = 0, the dynamics becomes that of 
parallel branes, with the action for the tachyon becoming just 



At late time, T is sufficiently large that we may neglect V{T). The Hamiltonian 




(44) 
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Figure 7: The momentum density Tqi as a function of t and x. The upward spike 
corresponds to the worldhne of the zipper with negative velocity while the trough 
corresponds to the zipper with positive velocity. 
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Figure 8: The pressure Tu as a function of t and x. Note that the pressure vanishes 
in the condensing region. 
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Figure 9: To determine the extent to which the late-time tachyon obeys (H6l) . we plot 
( !50l) as a function of x at fixed t = 6. This function approximately vanishes in the 
zipped region, showing that the dust description is valid. 



becomes just 

and the equations of motion ( l35i) and ( l36l) reduce to 



rp/Z 

UtT' 

/l + T' 



(45) 

(46) 
(47) 



This is the well-studied pressureless tachyon dust [24-28] (see also [29]) where the 
velocity of the dust is given by 

= -d,T (48) 



and the local rest energy density is 



(49) 



We find that the numerical solution accurately reproduces the tachyon dust once 
the branes have been zipped together. At x = where T' = 0, we saw in Fig. [3] that 
T — > 1 in agreement with ( l46l) . More generally, we plot in Fig. [9] the function 



1 - + T" 

1 + r2 + T'^ 



(50) 



which provides a normalized measure of the degree to which ( 1461) is satisfied numeri- 
cally. To good accuracy, fl3U|) is zero in the condensing region around x = 0, implying 
that in fact fHUj) is approximately obeyed. Furthermore, we see from Fig. [H] that the 
condensing region behind the zipper is almost exactly pressureless. 

If this zipped region is essentially a parallel Dl-Dl, we arrive at the following 
interpretation for the final configuration of the system. At late times there are two 
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curved Dl-branes separating from each other and straightening but connected to each 
other by the parallel Dl-Dl. Based on the results from the D-D description, the two 
separating Dl-branes is expected. However, rather than being devoid of branes, the 
region between them contains this decaying Dl-Dl pair. In the next section we will 
discuss the relation between these two descriptions. 

6 Relationship between D-D and D-D systems 

We have been analyzing the system of intersecting Dl-branes using a description in 
terms of a D-D pair. This choice of variables, at the level of an exact worldsheet 
description, amounts to an arbitrary gauge choice; a Dl is simply a Dl with opposite 
orientation, and a pair of Dl-branes intersecting with an angle 9 is equivalent to a 
Dl crossing a Dl at an angle if = ir — 6. However, the two effective field theory 
descriptions of these systems are inequivalent; each integrates out a different set of 
massive modes, resulting in different tachyonic and low-mass modes remaining. For 
example, using the YM approximation to the non-abelian DBI for the D-D pair, the 
tachyon's mass (I2l) matches the worldsheet value dH) only for 6* ~ 0, while in the 
tachyon DBI description of the D-D system, the tachyon's mass ( l30i) matches for 

6' ~ TT. 

Even though the two effective descriptions are not identical and the results need 
not to match a priori, as we will argue, they actually do produce dynamics which are 
the same, at least qualitatively. However, because the change of variables relating 
the two descriptions is nontrivially complicated, seeing that the results are physically 
similar can be a bit subtle. Before considering the relation between the D-D and 
D-D systems, we will begin with a simple illustrative example, where the 

two effective descriptions are in fact exactly equivalent and the change of variables 
between them is explicit. 

In the intersecting D-D system, how one describes the decay process depends on 
one's choice of gauge. In addition to choosing brane or anti-brane labels, the way 
the ends of the branes are connected is also to some degree gauge-dependent. The 
multiple equivalent descriptions can be easily seen from the non-abelian YM analysis 
of a Dl-Dl decay [15], as illustrated in Fig. [101 

In the initial variables, the matrix y giving the brane position is diagonal, corre- 
sponding to two straight Dl-branes. But, as the decay proceeds, off-diagonal elements 
begin to grow, fuzzing out the intersection point. However, y can be rediagonalized 
by a gauge transformation U which exchanges the brane ends [15,17]. In the new 
variables, the branes are disconnected, curved, and separating with time. Evolving 
back in time to the beginning of the decay, these two bent branes touch at a point. 
This odd configuration is, of course, physically inequivalent to the initial straight 
branes, because at t = the gauge transformation U relating them is singular. 

The important feature, which also appears in the D-D picture, is that in the 
original variables the decaying branes are not seen to annihilate at the intersection 
point, and the ends remain connected as in the initial state. Instead, the decay region 
is characterized by an unusual, non-intuitive state such as the stringy fuzz, and only 
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Figure 10: The process of D-D decay and recombination is shown in two gauges. On 
the left, the intersection of two straight branes develops off-diagonal elements and 
becomes fuzzy. After a gauge transformation U, the branes have reconnected and 
are pulling apart. However, because U is singular at the beginning of the decay, the 
initial conditions are inequivalent. 

after a gauge transformation does the system appear as two disconnected branes. 

For the description in terms of an intersecting D-D pair, we showed in Section 
[5] that the branes do not annihilate or reconnect into a disconnected pair of Dl- 
branes. Such an event is, in fact, outside the range of the tachyon DBI effective 
theory we used. Instead, the branes became zipped together, and where the Dl 
and Dl were parallel, the rolling tachyon was identified as inhomogeneous tachyon 
matter. However, although we do not have an explicit gauge transformation analogous 
to U in the D-D case, there is at least in principle a change of variables such that 
the parallel Dl-Dl are replaced by a gap between a pair of disconnected Dl-branes. 
Furthermore, in a more complete theory the parallel branes with tachyon matter could 
be alternatively described by a gas of closed strings and DO-branes (which themselves 
would decay) between the disconnected Dl-branes. 

Although we can not perform the explicit change of variables for the intersecting 
D-D pair, we can illustrate how it works in a similar but physically distinct system. 
For the initially straight intersecting Dl-Dl considered in Section [5l the relative 
position of the branes y was always antisymmetric, y{x) = —y{—x). We can just 
as easily solve the equations of motion but with y{x) = y{—x) instead, although for 
numerical computation we have to smooth out the branes at x = so as to avoid 
singular first derivatives there. This symmetric boundary condition corresponds to 
two angled branes with opposite orientations which touch at the point where they are 
bent. In addition, we will choose the intersection angle to be = 7r/2. This system is 
unstable, and there are two different directions in which it can decay corresponding 
to the two ways the branes can recombine, as shown in Fig. [TTlFI 

■^For the numerical computations, which decay mode the system takes depends on how exactly 
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Figure 11: The decay of two angled Dl-branes which touch at a point, such that y{x) 
is symmetric, can proceed in two ways. Final state (a) consists of a disconnected Dl 
and Dl. Final state (b) also contains a separating Dl and Dl, but they are connected 
by a parallel Dl-Dl. 



Figure 12: y as a function of x, for t = (black), 2 (dotted red), 4 (dot-dashed blue), 
and 6 (dashed green). 

One decay channel, mode (a) in Fig. [TTl closely resembles the late-time state of 
intersecting Dl-branes, except that instead of ending up with two separating Dl- 
branes, the orientation of one of the Dl-branes reversed. The numerical solution, 
presented in Fig. [T21 shows the Dl and Dl growing further apart while the tachyon 
field stays very close to zero. 

The other decay mode corresponds to the other way of recombining the ends of 
the branes and is (b) in Fig. [HI In the final state there is still one curved Dl and 
one curved Dl, but now they are connected. Numerically solving the equations of 
motion yields an evolution qualitatively very similar to those of Section the branes 
zip together and then the tachyon rolls towards infinity on the parallel Dl-Dl pair 
connecting the two separating Dl- and Dl-branes. 

Initially, these two final states appear quite different, but in fact, they must be 

the corner at x = is resolved. The details, however, are not that relevant to our discussion. 
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equivalent due to the symmetry of the initial state which is invariant under the 
exchange of x and y% This operation should map the final states into each other, and 
it does up to the parallel Dl-Dl. However, if we also exchange the roles of x and y 
in the tachyon DBI action, writing x and T functions of y, the parallel Dl-Dl would 
be present in decay (a) and not in mode (b). 

We conclude that the two modes are, up to this exchange of directions, different 
descriptions of the same physical state. In this example the change of variables is 
particularly simple; whether a parallel Dl-Dl connects the separating Dl and Dl or 
not depends on the choice of the worldvolume coordinates. At least at the level of 
open strings where we are working, both the parallel Dl-Dl with a rolling tachyon and 
empty space are equivalent. In particular, since at late time the parallel Dl-Dl with 
rolling tachyon does not support open string modes, it is effectively as if there are no 
branes there at all. 

Although the change of variables involved in the intersecting Dl-Dl system is 
much more complicated, the intuition gained from this symmetric example should 
hold. Because they are different approximations of the same system, the non-Abelian 
DBI D-D effective theory and the tachyon DBI D-D effective theory do not give 
identical quantitative results. However, since we should regard the parallel Dl-Dl as 
essentially the vacuum, both the D-D and D-D are qualitatively the same. 

7 Summary and discussion 

In this paper we studied the Dl-Dl configuration, initially intersecting at an angle 
ip. This configuration is unstable because tachyonic modes are present at all angles 
(except for the parallel D-D-brane pair (^9 = vr). We modeled the evolution of this 
system by deriving the equations of motion for the tachyon T(t, x) and the separation 
field y{t,x) from the tachyon DBI action and solved them numerically. 

We found that, at the very beginning of the evolution, the tachyon, which was 
localized at the intersection point, slowly rolled away from the maximum point of 
the effective potential. Then, after some time value, which could be interpreted 
as the (local) lifetime of the brane system [30] , the tachyon began to grow linearly 
with time and thus induced the dynamics for the separation field y{t,x). As time 
passed, the branes were pulled toward each other such that the point where they first 
met moved from to larger ±x values. The process resembled that of two zippers 
moving at opposite directions with speeds of light, zipping the branes together, and 
continuing indefinitely. The region between the zippers had the behavior of a decaying 
parallel Dl-Dl pair, but an inhomogeneous one. At late time, we were able to capture 
the essential features of the remnant by the well-studied inhomogeneous pressureless 
tachyon dust. 

Although the tachyon DBI is by construction accurate for ^ ~ tt while the YM 
description is valid for ^ ~ 0, both these results in terms of a Dl-Dl pair and the 
description in terms of intersecting Dl-branes give qualitatively the same evolution. 
The final state in both pictures contains two reconnected separating Dl-branes. We 

^The exchange corresponds to a 7r/2 rotation and an orientation reversal. 
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have argued that the parallel Dl-Dl connecting the branes in the tachyon DBI picture 
is an artifact of the choice of variables and is, in fact, equivalent to the gap between 
the branes in the YM description. 

In both of these effective theories the trace prescriptions for the actions are some- 
what ambiguous. In the case of D-D, the proposed action of [31] gives an alternative 
way of performing the trace and potentially gives different results from those found 
here. The non-Abelian DBI effective description of the D-D system should be valid 
not just at small angles. If the correct trace prescription could be found, for example, 
the worldsheet formula ([1]) for the mass at all angles may be reproduced. 

Another potential extension of these results is to other systems featuring localized 
tachyon condensation. One may have hoped this simple flat-space system could serve 
as a toy model for inhomogeneous brane decays in non-trivial curved backgrounds. 
In the Sakai-Sugimoto model, for example, a parallel D8-D8 pair decays into a single 
U-shaped D8 via the condensation of a tachyon localized at small radii [4]. However, 
one important qualitative difference between this and intersecting branes is that in 
the fiat-space case the decay process continues without end and the condensing region 
does not stay localized but instead grows without bound. 
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